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ABSTRACT 



We have made an extensive survey of emission-line stars in the IC 1396 H 
II region to investigate the low-mass population of pre-main sequence (PMS) 
stars. A total of 639 Ha emission-line stars were detected in an area of 4.2 
deg 2 and their i'-photometry was measured. Their spatial distribution exhibits 
several aggregates near the elephant trunk globule (Rim A) and bright-rimmed 
clouds at the edge of the H II region (Rim B and SFO 37, 38, 39, 41), and near 
HD 206267, which is the main exciting star of the H II region. Based on the 
extinction estimated from the near-infrared (NIR) color-color diagram, we have 
selected pre-main sequence star candidates associated with IC 1396. The age and 
mass were derived from the extinction corrected color-magnitude diagram and 
theoretical pre-main sequence tracks. Most of our PMS candidates have ages of 
< 3 Myr and masses of 0.2-0.6 M . Although it appears that only a few stars 
were formed in the last 1 Myr in the east region of the exciting star, the age 
difference among subregions in our surveyed area is not clear from the statistical 
test. Our results may suggest that massive stars were born after the continuous 
formation of low-mass stars for 10 Myr. The birth of the exciting star could be 
the late stage of slow but contiguous star formation in the natal molecular cloud. 
It may have triggered to form many low-mass stars at the dense inhomogeneity 
in and around the H II region by a radiation-driven implosion. 

Subject headings: stars: formation - stars: pre-main sequence - H n regions - ISM: 
individual (IC 1396) - open clusters and associations: individual (Trumpler 37) 
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Introduction 



As observational instruments have developed, systematic studies of stellar clusters 
have revealed many low-mass populations of pre-main sequence (PMS) stars in massive 
star-forming regions. It has been found that most stars are bo rn in groups or clusters, not 



in isolated regions (e.g., 



Lada fc Lada 



2003 



Allen et al. 



20071 1 . Stellar clusters provide a 



tool for investigating stellar evolution, because they are considered to constitute a natural 
sample of coeval stars at the same distance and with the same chemical composition. Young 

ution over a wide range of masses. However, 



19971 ) has revealed that the 



clusters appear to be a snapshot of stellar evo 
a careful study of the Orion Nebula Cluster ([Hillenbrand 
younger stars ( < 0.3 Myr ) are clustered towards the center of the Trapezium Cluster (a 
projected distance of < 0.3 pc from 8 1 C Ori), while the older stars are widely distributed. 
This indicates that star formation has recently occurred in the center of the cluster and 
that st ar formation with a low rate has taken place over a long time scale over the whole 



region. 



Muench et al. 



(120081 ) has suggested that there is no physical reason to separate 



the Orion Nebula Cluster and the Trapezium Cluster, and the entire region is a single 



contiguous star forming event. There are also examples of stellar c 



mass ive stars. Cep A HW2 is the brightest radio source in Cep A (IHughes fc Wouterloot 



usters still producing 



19841 ). an active star forming regio n within the molec u 



with an age of 5.5 Myr (Cep OB3b; 



Jordi et al. 



19961 ). 



ar cloud associated with the cluster 



Jimenez-Serra et al. 



(I2007h found a 



hot core internally heated by a massive protostar, and a circumstellar rotating disk around 
HW2. This suggests that a massive star is still forming in the young cluster. Furthermore, 
X-ray stars with a mid-infrared (MIR) spect ral energy distr i bution of young stellar objects 



(YSOs) around HW2 have been reported by 



Pravdo et al. 



20091 ). It is most likely that 



a low-mass PMS population was already formed before the formation of high-mass stars. 
Once a high-mass star starts to radiate significant amounts of UV radiation, the surrounding 
material is affected dramatically and the natal molecular cloud is disrupted and the birth 
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of the next generation of stars might be triggered. 



The apparent age spread o bserved in the color - magn itude diagram of young stellar 



groups has been investigated by 



Burningham et al. 



(120051 ). They found that the spread 



could not be explained adequately by the combination of the binarity, photometric 
uncertainty and variability of the members. In the case of the a Ori young stellar group, 

e explanation is t hat accretion from 



19991 ). However, another 



they suggested an age spread of 2-4 Myr. A possib 
a disk causes scatter of the stellar continuum light (ITout et al. 
possibility is that stars in a cluster are formed over a period of time. Three-dimensional 
numerical simulations including protostellar outflow feedback show that the cluster-forming 



clumps can maintain a quasi-virial equilibrium state for several 



rtNakamura fc Li 



periods of free-fall time 



Wang et al. 



L0) 



(2010) also 



20071 ). suggesting continuous star formation, 
demonstrated outflow-regulated massive star formation. The episodic star formation causes 
an intrinsic spread in age and hence in the position in the color-magnitude diagram. 



We have reported that young stars near the exciting star of the W5E H H 



region are 



20081). We 



systematically older than those near the edge of the H II region (INakano et al. 
further suggested that the birth of the most massive star occurs in the late stage of star 
formation in the cluster. Although the structure of the H II region is simple to investigate, 
the distance of W5 is large, at 2.3 kpc, and therefore it is not easily accessible to the wide 
range of stellar masses. 

IC 1396 is a ring-shaped giant H II region in Cep OB2 with a size of about 3°, located at 
the edge of the Ce pheus bubble, which is a large, far-infrared shell with a diameter of ~10° 



( iKun et al. 



20081 ) . The structure of IC 1396 is rather simple and is powered by the multiple 
O star system HD 206267 (RA = 21 h 38 m 57 ! :62, Dec = +57 o 29'20'. / 5), which is located near 
the geometrical center of IC 1396 in the young cluster Trum pler 37. The bright members 



of Trumpler 37 have been investigated by many authors (e.g., 



Garrison &: Kormendy 



1976 
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of bright-rimmed globules (IPottasch 



Marschall et al. I Il990t Ide Zeeuw et al.lll999t IContreras et al.l 120021 ) . There are a number 



19561 ) at its periphery. Its modest distance (< 1 kpc) 



and low foreground extinction (Av < 1 mag) make IC 1396 a good subject for our study of 
low-mass population associated with OB stars. The mean distance of 76 members of Cep 



1999|) 



Contreras et al 



0132 derived from the Hipparcos parallax is 615 pc (Ide Zeeuw et al 
(120021 ) made a color-magnitude diagram of the intermediate-mass members and fitted them 
to main-sequence stars at 870 p c. In this paper, we adopt 870 pc, which is frequently 



quoted by the recent works (e.g. 



Sicilia-Aguilar et al. 



20041 ) 



Ha emission, a characteristic of T Tauri stars (TTS), is radiated from the solar-type 
magnetic activity in weak-lined T Tauri Stars (WTTS) or the chro mospheric activity a nd 
magnetospheric accretion shock in classical T Tauri stars (CTTS) (jBriceno et al.l 120071 ). 
Thus, a survey of Ha emis sion-line stars is a po werful method for searching for young 



low- mass populations (e.g., 



Nakano et al. 



19951 ). Recently, NIR observations with array 



detectors and MIR observations from space (e.g. with the Spitzer Space Telescope) have 
advanced the protostar survey in the dense parts of giant molecular clouds. However, 
the Ha survey is still useful in areas with moderate or low extinction, and it is also a 
complementary method to search for PMS populations in a wide and rather dispersed 
region. Here, we report an extensive survey of emission-line stars in the IC 1396 H II region 
and discuss the formation sequence of low and high mass populations in IC 1396. 



Observations 



Slitless grism spectroscopic searches for Ha emission objects in IC 1396 were conducted 
in 2004 and between 2006 and 2009 using the Wide Field Grism Spectrograph 2 (WFGS2; 



Uehara et al. 



20041 ) at the f/10 Cassegrain focus of the University of Hawaii (UH) 2.2-m 



telescope on Mauna Kea. The log of the observations is given in Table 1. A 300-line 
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mm 1 grism was used in conjunction with a wide Ha filter to isolate the first-order spectra 
between 6265-6765 A. A Tektronix 2048x2048 CCD was used as the detector, yielding a 
dispersion of 3.8 A pixel -1 , a pixel scale of 0'.'34 and an instrumental field of view of 11(5 x 
11(5 

The whole IC 1396 region was covered by 157 field positions including overlapping 
areas. About 4.2 deg 2 were observed, both in the grism and the direct imaging mode of 
WFGS2. Figure 1 shows the boundary of the observed fields and indicates the names of 
bright-rims, dark clouds, and two O stars in the IC 1396 H II region. For each field, we have 
taken a single 30-s exposure in the direct mode with a wide Ha filter. Then, we inserted 
a grism and obtained three dithered frames of 300-s exposures for slitless spectroscopy. 
Direct images with an i'-band filter were also taken with WFGS2 in 2008 August and 2009 
August. Three 30-s exposures were made for each field. The i'-band flux is less affected 
by excess emission of the circumstellar disk, Ha emission, and continuum veiling of CTTS. 
Thus, z'-band photomet ry is useful as a measure of the underlying photosphere of CTTS 



Cieza et 



( jSmith et al. 



al. 



20051 ). 3-5 standards for the Sloan Digital Sky Survey (SDSS) system 



20021 ) were observed each night at nearly the same airmass (difference of < 



0.1) as the targets. We used a dome flat for flat fielding. 

The data were reduced using the IRAF software package. We applied the standard 
procedure for bias subtraction and flat fielding. For each field, three dithered grism frames 
are combined into one frame by referring the position of the Ha emission or absorption 
line of the star(s) in the field. We picked all stars in each wide Ha image using the 
DAOFIND package and extracted their one- dimensional spectra from the grism image. We 
visually investigated the spectra for the presence of the Ha emission line and found 639 Ha 
emission-line objects in the whole surveyed area. We used the SPLOT package to measure 
the equivalent width of the Ha line in the spectra. Their coordinates were determined 
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by the astrometric calibration tool of the Starlink program GAIA using USNO-B catalog 
stars as the reference in the field. All of the photometric measurements were made with 
the APPHOT package for aperture photometry. The typical seeing of the image was 
2.7-3.0 pixels. We used a 4-pixel diameter aperture for the aperture photometry. We 
then calibrated the photometric values to the SDSS system via standards. In general, the 
run was of high quality, but thin cloud affected a few fields. We checked the stability 
of the magnitude for three direct images in each field and examined the quality of the 
measurement. Excluding stars with saturated images or images with low photometric 
quality, we measured 586 stars in 13.6 < i' < 19.3. 

3. Results 

3.1. Spatial Distribution 

We list the 639 Ha stars in Table 2. The table shows the (1) identification number, 
(2, 3) J2000 coordinates, (4) z'-band magnitudes and (5) equivalent widths of Ha emission 
(A). The Ha equivalent widths were measured for 608 stars. We compared our list of 
emission stars with the Two Micron All Sky Survey (2MASS) Point Source Catalog with 
12 stars supplemented by the 2MASS 6x Point Source Catalog. We found matches for 
617 stars with 2MASS stars within 2" and determined that 548 of these have good quality 
JHKs photometry (quality flag AAA). Due to small separation angles or positional 
uncertainties, two pairs of Ha stars are identified with single 2MASS sources (No. 250/251 
with J21381703+5739265, and No. 600/601 with J21451578+5717392). We have included 
these stars in the above number. We list the 2MASS designation of the NIR counterpart, 
photometric quality, J magnitude, and J-H, H-Ks colors in columns (6)-(10) of Table 
2. Av estimated from the NIR color-color diagram (see §4.1) is also listed in column 
(11). About 70% of our Ha stars cannot be identified with known sources in the SIMBAD 
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database. The identification in the literature or comments are listed in column (12) 



Figure 2a shows the spatial distribution of emission-line objects while Figure 2b 
represents the contour map of the number d ensity of the detected emission-line objects 



with the kernel method ( Gomez et al. 



1993|). We adopted a Gaussian shape for the density 
distribution of the kernel with a smoothing parameter of h=2'4. The isolines are drawn 
at intervals of 150 stars per deg 2 , which is about three times the average number density 
in the field region (see Table 4), starting from 100 stars per deg 2 . It is clear from Figure 
2 that the distribution of emission-line stars is far from uniform. To demonstrate the 
different properties in the different areas, we divided the surveyed field into the following 
subregions according to the number density of Ha stars. The primary peak of the most rich 
aggregate (the West region) is located between the main exciting star HD 206267 and the 
western elephant trunk globule Rim A (IC 1396 A). The secondary peak is at the opposite 



side, to the east of HD 206267 (the East 
the bright-rimmed globule SFO 38 (SFO: 



region), and the thir d peak is at the south of 



Sugitani et al. 



1991 



the SF038 region). Other 



small aggregates of emission-line stars are located near SFO 37 (the SF037 region), SFO 39 



and SFO 41 (the S 
365 dark cloud (B: 



'039/41 



region), Rim B (the Rim B region), and at the north of the B 



Barnard 



19271 ) or around (21 35 m , +57°). See Figure 1 for references. 



The West and the East regions appear to be surrounded by a halo-like structure (the Halo 
region). Some Ha stars appear to be associated with other small bright rimmed clouds such 
as Rim C and SFO 34 (Rim D). We defined the rest of the region as the field area (the 
Field region). 



3.2. Source Identification 



Our sample of Ha stars contains two carbon stars CGCS 5454 (No. 575) and CGCS 
5401 (No. 287), and one Be/X-ray binary Cep X-4 (No. 340), whose distance was confirmed 
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to be 3.8 kpc by iBonnet-Bidaud fc Mouchet I (119981 ). Our sample of Ha stars also contains 
a candidate of the exciting so urce of the Herbig -Haro object HH 588 (No. 436), which is 



deeply embedd ed in SFO 37 (llkeda et al. 



426 (No. 236) flSicilia-Aguilar et al. 



trans i ent objects ( 



2004 



Sicilia-Aguilar et al. 



Mercer et al. 



2009 



20081 ). one Herbig Ae/Be star (HAeBe) MVA 



2007), an d ~ 100 TTS including Class-II, III, and 



2006b 



Barentsen et al. 



2007 



201lh . 



Getman et al. 



2007 



Reach et al. 



An objective prism survey using the Schmidt telescope 



Kun 



(|1986h detected 155 



Ha emission objec ts in IC 1396 and an additional 65 emission stars were reported in a 



subsequent paper ( IKun fc Pasztor 



19901 ) . Their spatial distribution shows no noticeable 
concentration or clustering. Out of 220 Kun Ha stars, 150 stars are located in our observed 
field; however, only 11 stars (No. 287 = Kun 80, No. 13 = Kun 171, No. 340 = Kun 193, 
No. 618 = Kun 215, No. 49 = Kun 308, No. 80 = Kun 310, No. 97 = Kun 314 S, No. 252 
= Kun 321, No. 318 = Kun 325, No. 410 = Kun 330, No. 487 = Kun 332) are identified 
among our Ha stars. Although Kun 315 (LkHa 349) is not included in Table 2, our grism 
spectrum appears saturated but indicates a P Cyg-type Ha profile. Kun 320 (GL Cep) 
and Kun 337 (HD 239745) are also too bright and saturated in our grism spectra, and the 
positions of Nos. 77, 230, 12, and 266 are offset from Kun 53, 73, 170, and 322 by around 
10". Even if we include these seven stars, only 12% of stars can be identified with Kun 
stars. The rest of the Kun stars show obvious Ha absorption or no emission, and 10 of them 
show the spectra o f M-type stars in our grism spectra. Medium resolution spectroscopy by 



Balazs et al. 



(Il996[ ) show that Kun 193 (No. 340; CepX-4) exhibits conspicuous hydrogen 
emission, but 34 other Kun stars d o not. They suggested that the lack of emission is due 



to the time variability of Ha stars. 



Contreras et al. 



( 120021 ) failed to detect Ha emission in 



spectroscopic observations of 20 selected Kun stars, only finding emission in Kun 314 S. 
Thus, they suggest that the original identification of emission was incorrect. Our result is 
consistent with their suggestion. 
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Sicilia-Aguilar et al. I ( 12004 . 120051 ) extensively studied the low-mass population of the 
two clusters in Cep OB2, Trumpler 37 and NGC 7160, to investigate the disk evolution at 
ages 1-10 Myr. For Trumpler 37, they have performed optical photometry and multifiber 
spectroscopy in a square about 45' per side, centered at HD 206267. Approximately 45% of 
the low- mass population indicate MR exce sses, suggesting heated dust in the circumstellar 
disk. Finally, ISicilia-Aguilar et al. I (j2006b| ) examined the membership using Ha emission, 
Li A6707 absorption lines, and radial velocities as membership criteria. They identified 
accretion properties in 170 members in the central part of IC 1396. 53% of the members of 
this list are also detected in our survey. The equivalent width of th e Ha emission-line of 



our candidates is consistent with that of 



Sicilia-Aguilar et al. 



( 120051 ). except for only a few 



stars. We have detected ~ 90% and < 30% of Class-II and Class-Ill objects in their list as 
Ha stars. This result is naturally explained by the higher detectability of weak emissions in 
their high dispersion slit spectroscopy. 



Using archival data of Chandra and Spitzer, iMercer et al. I (120091 ) have identified 25 
X-ray sources as cluster members, which have MIR colors of CTTS or WTTS in the 10' x 
8' region around HD 206267. Their sources consist of five CTTS and 20 WTTS. Although 
four of the CTTS were detected in our observations, only two WTTS match our Ha stars. 



Getman et al. 



(120071 ) obtained X-ray images of the bright rimmed cloud SFO 38 (IC 
1396N) using the Chandra ACIS detector. They identified X-ray sources in the area of 
a square of 5' with the Spitzer IRAC sources and concentrated their study on 25 sources 
associated with the globule SFO 38. Of the 117 X-ray sources, 17 sources are identified 
with our Ha stars and five of these X-ray sources are associated with Ha stars in the SFO 



38 cloud. 
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3.3. Comparison with IPHAS Sources 



In the recent paper by lBarentsen et al. I (120111 ). they identified 158 PMS candidates in a 
7 deg 2 area towards IC 1396 by using the data of IPHAS (t he Isaac Newton Telescope/ Wide 



Field Camer a Photometric Ha Surve y 
identified by 



Drew et al. 



Barentsen et al. 



2005|). Out of 143 PMS candidates 



(120111 ) in our field of view, 119 stars are in common with 
our Ha sample. We examined the common Ha stars in IPHAS to confirm our photometric 
quality. The IPHA S magnitude is based on the V e ga magnitude and henc e we converted it 



to AB magnitude (IGonzalez-Solares et al. 



20081 ) 



Barentsen et al. 



( 1201 ll ) restricted their 



sample to sources with r' IPHAS < 20 mag, which corresponds to i' < 19 mag if we adopt 
the colors of known T Tauri stars in their figure 4. For most of the stars in common, our 
^'-magnitudes are consistent with those in the IPHAS catalog within 0.1 mag. Although 
a large scatter (a = 0.3 mag) of the magnitude differences is found, it may be due to the 
photometric variability of young stars. 



Barentsen et al. 



Sicilia-Aguilar et al. 



2011 



compared their objects with the members confirmed by 



(j2006b[ ). and esti mate the completeness . 89% o f PMS candidates 



with the Ha EW > 30 A identified by ISicilia-Aguilar et al. 



Barentsen et al. 



(l2006bf ) are recovered by 



(120 111 ). We also estimated the completeness of our survey with the same 



manner. F or the PMS candidates in cluding; smaller EW ( > 10 A ), the fraction of recovered 



sources by 



Barentsen et al. 



(120111 ) falls off to 50%. On the other hand, the fraction of 



recovered sources with EW > 10 A maintains 86% (51 out of 59) in our sample . We 



Barentsen et al. 



(2011) 



confirmed that none of their rejected candidates ( 30 stars ) by 
were included in our list. Thus, our observations are more sensitive than their survey for 
the emission-line stars with small Ha equivalent width. 
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3.4. Identification with AKARI Sources 



We identified 27 sources among our Ha stars with the sources of the point source 



catalog of the MIR A 



1-Sky Survey obtaine d with the infrared camera (IRC) on board 



the AKARI satellite ( llshihara et al. 1120101 ). with position differences of less than 5'' 
Table 3 gives the AKARI 9-/zm (S9W) and 18- /jm (L18W) magnitudes and their 2MASS 
photometric data. Their positions on (J-Ks) vs. (S9W-L1SW), (J-Ks) vs. (Ks-S9W), 



and (J-Ks) vs. (Ks-L18W) diagrams are shown in 



TTS from these diagrams, proposed by 



Takita et al. 



igure 



3. The criteria for extracting 



( 20101 ). are also shown dashed 



line (Figure 3a) or areas (Figure 3b, c). The colors of three stars (Nos. 287, 309, 575) are 
consistent with cool carbon stars, and about 20 stars are pre-main sequence stars. The 
bright source No. 436 is the candidate for the exciting source of the HH 588 object. Figure 
4 shows the contour map of the number of emission-line stars overlaid on a false-color map 
using the IRC on AKARI. Blue is for 9 //m and red is for 18 /im. Probable members of 66 



■ 3-F stars from IContreras et al 



deZeeuwetal 



(120021 ). supplemented by 8 OB stars from Hipparcos results 



19991 ). are also shown as filled blue circles. As for Ha stars in the West 



and East regions, high or intermediate-mass stars are concentrated near the center of the 
large dust cavity, which is presumably swept up by the expansion of the H II region. We 
note that their distribution is slightly shifted to the east of HD 206267. 



4. Analysis & Discussion 

4.1. Probable Members 

Our list of Ha stars in Table 2 may contain foreground or background, magnetically 
active field stars. We estimated the number of unrelated main sequence stars with enhanced 
chromospheric and coronal activity. Within the 4.2 deg 2 field towards the IC 1396, the 
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Besancon model of the Galactic stellar populations (IRobin et al. 1120031 ) predicts ~4000 and 
~ 1700 dM stars down to V = 20 mag in the foreground and the background (the distance 
range between 870 and 1500 pc) field, respectively. Although 44% of the foreground dM 
stars are in the range of V < 19 mag, 3% of the background dM s tars a re in the same 
magnitude range. Following the same approach a s IDahm fc Simonl (120051) . we used the the 



stellar luminosity function derived by Hipparcos ( jJahreiB fc Wielen 



1997) and the dMe 



Hawley et al. 



(11999 1 



(with an Ha EW > 1 A) incidence from the nearby star survey of 
We adopted a completeness limit for the present observations of V = 20 mag. We assumed 
the scale height of 325 pc for an exponential falloff in the space density, and 0.7 mag kpc -1 
for the diffuse absorption. The estimated number of dMe stars in the foreground and the 
background (up to 1500 pc) is 193 and 33 in our observed area, respectively. The total 
number of the field dMe stars, 54 per square degree, is consistent with the average number 
density of Ha stars in the field area, 58 per square degree. 

It is also possible that there is contamination by pre-main seq uence stars associated 



with the foreground cloud (jPatel et al. 



1995 



Weikard et al. 



19961 ). As the AKARI 9 ^m 



image, which traces the photodissociation region in the periphery of the H II region (Figure 
4), is likely to be dominated by polycyclic aromatic hydrocarbons (PAHs) emission features, 
Ha stars associated with 9-fim emission are not likely to be foreground objects. The dark 



clouds B 163 and B 163SW do not show an y sign of brig ! 



it rims and their radial velocities 



of the molecular emission line are 6 km s 1 (jPatel et al 



19951 ). The radial velocities of the 



other associated glo bules are between —9 km s 1 and +2 km s . From the radial velocity, 



Leung et al. 



(119821 ) suggested 150 pc for the di stance of B 161, whic h also shows no bright 



rims and no indic ation of heating 



dark lane Kh 161 (jKhavtassi 



19961 ). The broad, curved 



rom behind (IWeikard et al. 
I960) and inverse S-shaped B 365, both of which show weak 
and diffuse emission at the northern and southern part of Figure 1, appear in foreground 
clouds. 
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Figure 5 illustrates the NIR color-color diagram of our emission-line stars. The 



solid lines yield t 



Bessell fc Brett 



le lo cus of the main-sequence and of the giant branch according to 



(119881 ) . whil e the dash-do t ted st raight line indicates the locus of the 



unreddened CTTS follow ing 



equations of 



Cutri et al. 



Mever et al. 



reddening vector adopted from 



(119971 ). We used the color transformation 



(I2008h to convert t 



neir 



Cohen et al. 



oci into the 2MASS system. The 



(I198ll ) is also shown. We selected nearly 40% 



of the Ha stars whose NIR colors are consistent with the colors of the reddened CTTS in 
Figure 5 and calculated the source extinction by tracing back along the reddening line to 
the CTTS locus. In the above procedure, we have not separated between WTTS and CTTS 
candidates, and treated all Ha stars as the same way. Although generally an unreddened 
WTTS does not show the color of the CTTS but rather of the normal main sequence, most 
WTTS have NIR colors of K- to early M-type stars, i.e., (H-Ks, J-H) ~ (0.2, 0.64), that 
is close to the locus of unreddened CTTSs. Thus, we have not lost many WTTS in the 
above procedure. Three carbon stars, discussed in the previous section and shown as open 
circles in Figure 5, are excluded. Figure 6 shows a histogram of the measured extinction of 
252 stars excluding three carbon star candidates (see §3.3) and No. 250/251. The median 
value of the source extinction is Av= 1.74 mag. This value is con sistent with the average 



extinction Av=1.56 mag obtained by 



Sicilia-Aguilar et al. 



( 120051 ) for Trumpler 37. 



The high mass stars in our surveyed area are t wo O-type stars, HP 2 06267 (06) and 



HD 206183 (09.5V), and six early B stars (B0-B5) f ide Zeeuw et al 



the fo r eground extincti o n of t hese OB stars should be 1.3-2.1 mag ( jGarrison fc Kormendy 



1999) . When Rv=3.1, 



19761 ). 



Contreras et al. 



(120021 ) excluded stars with Av < 1.0 mag as foreground stars and 



excluded stars with Av > 2.3 mag as background objects. As the background contamination 
by field dMe stars is not so high, many of the emission-line stars with Av > 2.3 mag in our 
survey appear to be associated with dark clouds or globules. Therefore, we selected 189 
stars with Av > 1.0 mag as likely members associated with IC 1396. Most known T Tauri 
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stars (e.g., ISicilia-Aguilar et al. Il2006bl ) in Table 2 have Av > 1.0 mag, but 12 T Tauri 
stars have low extinction, and more than half of them are located near HD 206267. 

Figure 7 shows a histogram of the measured equivalent width. The fraction of weak 
(EW < 10 A) emission-line stars among the IC 1396 PMS candidates is low (27%) relative 
to that over the entire sample (44%). It appears that the fraction of CTTS candidates 
among the PMS candidates in the IC 1396 region is higher than that among the candidates 
proposed as foreground (or background) objects. We could not find any evidence of the 
difference between the spatial distribution of CTTS and WTTS. 

To identify the HAeBe stars, which have larger infrared exces s than CTTS, we a dopted 



the extended criteria in the NIR color-color diagram proposed by iLee fc Lim I (120081 ). They 
used the region between two parallel reddening lines on their figure 4 and above (J-H) 
> 0.2 mag to pick up HAeBe stars, excluding B[e] and classical Be stars from the young 
star sample as contaminants. According to their criteria, we selected five stars, redder than 
the reddening line passing through the trun cated point of the CTTS locu s and J-H > 0.2 



mag, including a known HAeBe, MVA-426 (ISicilia-Aguilar et al. 



2006b|). It is suggested 



that from their J magnitudes Nos. 13, 236, and 606 are HAeBe stars and No. 479 and 545 
are Class-I objects. In the following discussion, we have added these stars and 12 T Tauri 
stars with low extinction to the probable members totaling 205 stars. The number of new 
unpublished PMS candidates identified in the current survey is 72 including three Class-I 
objects or HAeBe stars. 

Figure 8 illustrates the contour map of the number density of emission-line stars 
associated with IC 1396. Dashed boxes represent the eight subregions, namely, the East, 
West, SF038, SF037, SF039/41, Rim B, Halo, and Field regions. Their distribution is 
not very different from Figure 2, though the peak of the distribution at the East region is 
slightly shifted westward. While the aggregate at the north of B 365 (see Figure 1), which 
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is shown in Figure 2b, disappears in Figure 8 and therefore could be a foreground object, 
other aggregates should be associated with IC 1396. Table 4 summarizes the number of 
emission-line stars in each region. The Ha stars, whose Av values we could not obtain, 
appear to be spread over the observed field rather uniformly. Thus, we suggest a large 
number of them are field stars. 



4.2. Ages and Masses of Ha stars 



Marschall et al. I ( 119901 ) determined that the time after the fi r st sta rs reach the 



main-sequence is 7 Myr in Trumpler 37 and 



Sicilia-Aguilar et al. 



(120051 ) confirmed 4 



Myr for the age of Trumpler 37 from optical photometry and theoretical isochrones. The 
expansion of the H II region creates a swept-up shell a nd results in a comp ressed molecular 



Patel et al. 



( 119951 ) suggested a 



ring around its periphery. From the size of the ring, 
dynamical age of the system of globules of 2-3 Myr. 

Figure 9 shows the extinction-corrected i! vs. (i'—J) diagram of the 194 probable 



memb ers with z'-band data, using the i'-band extinction co rrection from 



( 119891 ). Here, we adopted the pre- main sequence model of 



Siess et al. 



Cardelli et al. 



tort . Their 



effective temperature and lu minosity are convert ed into our system using the table of 



Kenyon & Hartmann 



Jordi et al. 



(J2006J). In Figure 9, the 0.1, 



(119951 ) and the equation of 
1, 3, 5, 10, and 100 Myr isochrones are overlapped as well as the evolutionary tracks for 
masses from 0.1 to 2.0 M , at an assumed distance of 870 pc. As the reddening vector runs 
almost parallel to the isochrones, the extinction value would not severely affect our age 
estimate. 

The uncertainty of the distance modulus (8.94-9.70) and the choice of PMS models 
affects the estimation of ages and masses of Ha stars from Figure 9. We used the pre-main 



-re- 



sequence models of |P 
system by a table of 



Kenyon fc Hartmann 



alia fc Stahler I (|199S ) , whic h are converted to our photomet ric 



f|1995T ) and the eq 



for comparison. By assuming the nearest distance, 615 pc fide Zeeuw et al 



uation of 



Jordi et al. 



fl2006f ). 



1999|), the 



evolutionary tracks shift upwards, and it makes the age ~0.2 dex older in both models. 
The age and mass histograms of the probable members in different choices of PMS models 
adopting 870 pc are shown in Figure 10. They indicate that the numbers of young stars 
( < 1 Myr ) predicted from the Siess models (Fig 10a) are generally smaller than those 
predicted from the Palla & Stahler models (Fig 10b). On the other hand, the number 
of low mass stars ( < 0.4 M ) derived from the Siess models (Fig 10c) are slightly large 
compared to those derived from the Palla & Stahler models (Fig lOd). We should notice the 
incompleteness of ou r PMS candidate sam ple along with possible effects on the distribution 



of masses and ages. 



Da Rio et al. 



( 120101 ) carefully investigated the stellar population of 
the Orion Nebula Cluster by using these two PMS models. The completeness function in 
the mass-age plane, derived from their statistical simulation for two evolutionary models, 
shows that a large incompleteness correction is needed in the low mass (< 0.2 M ) range. 
They attributed the difference of results to the different shapes of tracks and isochrones 
in this mass range. Although we adopted the Siess models, the uncertainties would be 
larger for stars with ages < 1 Myr and masses < 0.2 M . Furthermore, we could not detect 
emission-line stars near the region of higher obscuration such as the globules, and within 
V of the bright O star HD 206267, due to saturation of the detectors, light contamination 
from this star, and diffraction spikes from this star. 

Table 5 shows the number of Ha stars in each age and mass bin in Figure 9. The age 
and mass histograms of Ha stars in each subregion are shown in Figure 11. We combined 
the SF037/38/39/41/Rim B regions in the BRC (bright-rimmed cloud) region. Most of 
the Ha stars have ages of < 3 Myr and masses of 0.2-0.6 M Q . The number of stars with 
an age less than 1 Myr and with an age of 1-3 Myr are comparable in the region associated 
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with bright rims, i.e., the BRC region. We confirmed that the above result is also true for 
the age assigned from the Palla and Stahler model, and is not dependent on the choice of 
distance. The histogram in the West region is similar to that of the BRC region, probably 
because most of the Ha stars in the West region are associated with Rim A. 



From the extensive survey of young stars by Spitzer, 



Koenig et al. 



(120081 ) found 



two distinct generations of star formation in W5. They consi dered two triggere d star 



formation mechanisms, radiative driven implosion (RDI; e.g. 



Miao et al. 



20091) and 



collect-and-collapse, to explai n their findings and conc luded that both mechanisms are at 



work in the W5 H II region. 



Deharveng et al. 



( 120 lOl ) extensively studied the nature of 



infrared bubbles which enclose H II regions and discussed star formation triggered by the 
expanding H II region at their edges. They show that most of these bubbles enclose H 
II regions ionized by 0-B2 stars, many are formed to be surrounded by cold dust shells 
(good candidates for C&C processes), and many have bright-rimmed dust condensations 
protruding inside the H II region (good candidates for RDI process). In the IC 1396 region, 
the presence of aggregates associated with the bright rims immediately inside the IC 1396 
H II region or the near side of the ionizing star HD 206267, such as the West (Rim A), 
Rim B, and SF038 regions, suggest that the star formation in these regions has advanced 
away from HD 206267 in < 3 Myr. Such a configuration is consistent with squeezing of the 



pre-existing condensa tions compressed by t 



structure of SFO 37 ( Su 



38 ( iGetman et al. 



gitani et al. 



1997 



re pressure o 



Ikeda et al. 



the ionized gas. The kinematical 



20081 ) and an X-ray study of SFO 



20071 ) also favored the RDI model. 



4.3. Star Formation in the Central Region 



In our previous results for the W5E H II region ( INakano et al. 1120081 ). we have found 
three aggregates of Ha emission-line stars in the W5E H II region. Two of them are located 
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at the edge of the H II region, and one is near the central exciting 07 star. The former two 
aggregates are systematically younger than the latter. 

Based on a comparison of the stellar ages in the BRC region with those in the East 
region of IC 1396 (Figure 11), it might appear that the low- mass population in the East 
region is older. We only find two stars younger than 1 Myr in the East region, in contrast 
to the West region. It may suggest that the star formation activity in the East region has 
ceased recently. As there remains no molecular material associated with the population in 
the East region, the exciting star HD 206267 has already dissipated the natal molecular 
cloud. However, a Kolmogorov-Smirnov (K-S) test gives a probability of > 15 % that the 
two samples have the same distribution in age, even if we choose either PMS model. The 
difference between the ages of Ha stars in the BRCs region and the East region is not 
statistically significant. 



Sicilia-Aguilar et al. I (l2006al ) estimated ages for six members near Rim A or in the 
globule. They found that the ages of five stars are younger than 1 Myr, which is significantly 
lower than the average age of Trumpler 37. Figure 12 plots the spatial distribution of 
Ha stars with derived ages. At a glance, younger stars in Figure 12a are preferentially 
located at the West region compared to those at the East region. Again we used K-S test 
to compare the distribution in right ascension and declination derived from different choice 
of PMS models. A probability that two samples of different age (Fig 12a and 12b) have the 
same distribution in two coordinate axis is > 25 % and > 54 %. Thus, the age difference 
of young stars in the West and the East region is not statistically clear from the spatial 
distribution of our sample. 

In addition to the West region, Figure 12b shows a number of stars with an intermediate 
age (1-3 Myr) in the East and SF038 regions, extending to the Halo region. A linear 
sequence of young stars across the center is noticeable. A filamentary dark cloud extends 
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from Rim A to the western rim of the H II region. It is noteworthy that Ha stars are 
distributed along this structure and further extend to the east of HD 206267. This spatial 
distribution strongly suggests a fossil record of a star forming cloud including Rim A and B. 



Sicilia-Aguilar et al. 



(120051 ) argued a spatial east-west asymmetr y in the young 



low-m ass members of Trumpler 37, compared with more massive ones (IContreras et al. 



20021 ). The same spatial asymmetry with respect to the central io nizing O star, HD 206 267, 



was se en in the distribution of the Spitzer IRAC excess objects (jSicilia-Aguilar et al 



2006al ). The asymmetry of Ha stars with high disk fraction may be explained by 



photo-evaporation of the di sk material around the forming stars by strong UV radiation. 



Sicilia-Aguilar et al. 



(120051 ) suggested that the interstellar material, lying between powerful 



06 stars and young stars with disks, shields the disk from UV radiation. How ever, there is 



no ev idence of such an asymmetry of interstellar ma terial from molecu 



19951 ) and millimeter dust continuum observations (jRosolowsky et al 



ar lin e (jPatel et al. 



201(1 ). Actually, 



significant destruction of disks occurs for stars located out to distances of < 0.5 pc, which 



corre s ponds to only a few arcmin 



20070 . 



Sicilia-Aguilar et al. 



rom HD 206267, as found in NGC 2244 (IBalog et al. 



(120051 ) additionally suggested probable age differences for the 
explanation for the asymmetry in the CTTS distribution. Instead of disk photo-evaporation, 
we prefer that this spatial asymmetry represents the star formation sequence. If most of the 
young low- mass population in the West region was formed as a result of a triggered event, 
those in the East region could be formed along with massive stars. 



4.4. Star Formation in the BRC Region 

There are 116 probable members of emission-line star in the BRC region including 
the West region, and the total mass is about 60-70 M . If we add 35 emission- line stars 
in the East region, the total mass reaches 70-90 M Q . Since our survey is biased toward 
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the region of low extinction, the above estimates of the number and mass of the low-mass 
population in the BRC region including the West and th e East region should be considered 



the lower limits. For example, 



Choudhury et al. 



( I2010l ) found ~45 YSOs in and around 



SFO 38, where we could not detect any Ha stars in the globule. Of these, 26 sources are 
Class-O/I, Class-I, and Class-I/II, but there are no such objects situated outside the ionized 
rim. It follows that the birth of low-mass stars associated with bright rims appears to be 
the primary mode of star formation in IC 1396 at present. 



The higher mass O-F stars of Trumpler 37 ( jde Zeeuw et al 



2002 


) total 


of 


Kroupa 



1999 



Contreras et al. 



2002 ) total ~ 200 M in our observed fields (Figure 4). Assuming the initial mass function 



(120011 ) between 0.01 M and 50 M , we estimate the total mass of stars in 
Trumpler 37 to be ~ 500 M Q . Young low-mass stars, which probably formed by a triggered 
mechanism in the BRC region, contribute a fraction of the mass of stars in IC 1396. 



4.5. The Birth of HD 206267 



Ikeda et al. 



(120081 ) has determined the alignment of YSOs and their age gradient at 
the tip of SFO 37. They argued that this alignment and the elongation of the head of the 
globule are attributed to UV radiation from the 09.5V star, HD206183, at their formation 
stage. The dominant exciting source, however, is currently HD 206267. Consequently, HD 
206267 formed at the late stage of st ar formation in the cent ral part of the H II region and 



disrupted the natal molecular cloud. 



Barentsen et al. 



(120111 ) reported the age gradient of 



the PMS candidates along the distance from this hot star. Their results also favors the 
triggered formation of stars in IC 1396. HD 206267 is known as the Trapezium system of 
OB stars with a maximum separation of 1 5000 AU, and the main component A is a triple 



(jStickland 



19951). 



Burkholder et al 



(119971 ) confirmed the observed masses and evolutionary 



model masses and suggested that the two bright components of the triple are 06.5 V((f)) 
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(Ai) with 38-61 M and 09.5: V (A 2 ) with 17-28 M . T herefore, the age of HP 2062 67 is 



estimated to be < 3 Myr from the main sequence lifetime (jSchaerer fc de Koter 



19970 . 



Numeri cal simulation consideri ng p rotostellar outflow-drive n turbulence has been 



discussed by 



Li fc Nakamura 



()2006j) and 



Nakamura fc Li 



(120071 ). Following their study, we 



suggested that the massive exciting star of the W5 E H II region was possibly formed after 



the continuous production of stars (INakano et al 



20081 ). The birth of the most massive 



star should be the last stage of star formation. The contiguous formation of stars for 10 
Myr in Trumpler 37 could be supported by Figure 12. Many of the low-mass stars near 
Rim A and SFO 38 might be formed in < 3 Myr. As there is no conspicuous molecular 
material except the elephant trunk globule, Rim A, in the evacuated dust cavity (Figure 4), 
HD 206267 was likely born in the central part after the continuous low mass star formation. 
It is reasonable to suppose that the formation of massive stars and the associated low-mass 
population in Trumpler 37, continued for a few Myr. Although it is difficult to separate the 
low-mass stars formed as the same generation with HD 206267 and formed by the trigger 
mechanism of HD 206267, we believe that many of H a stars at the rim were triggered to 
form by a radiation driven implosion mechanism after the formation of HD 206267. 



5. Summary 

We have made an extensive survey of emission-line stars in the IC 1396 H II region 
to investigate the low-mass PMS population near the OB stars. A total of 639 Ha stars 
were detected over a 4-deg 2 area and their i'-photometry was performed. Their spatial 
distribution shows several aggregates near the bright-rimmed clouds at the center and the 
edge of the H II region (Rim A/B and SFO 37, 38, 39, 41) and one near the main exciting 
star of the H II region. Based on the source extinction estimated from the NIR color-color 
diagram, we have selected PMS candidates associated with IC 1396. The age and mass were 
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derived from the extinction corrected color-magnitude diagram and theoretical evolutionary 
tracks. Most of the Ha stars have ages of < 3 Myr and masses of 0.2-0.6 M . 

The spatial distribution of the PMS candidates indicates that the stellar population 
in the East region is spatially distinct from that in the West region, associated with 
Rim A globule. The age dif f erence between two pop ulations was also suggested by 



Sicilia-Aguilar et al. 



f l2006bh : 



Barentsen et al. 



( l201lh . The spatio-temp oral gradients o 



the stel 



( 1200 



ar age were also d e monstrated near th e brig 



Ikeda et al. 



teOQSh : 



Choudhury et al. 



it-rimmed clouds by 



Getman et al. 



( 1201 Oh . Our K-S test does not support the 



significant difference between ages of two populations. Also, the spatial distribution of < 1 
Myr stars is not significantly different from the distribution of 1-3 Myr stars. Although the 
star formation activity in the East region could be low in the last 1 Myr, we need further 
confirmation. 

The linear sequence of stars including the West (Rim A) region and Rim B region, 
which is conspicuous in Fig 12b of 1-3 Myr age across the center of the H II region, suggests 
a fossil record of the star forming cloud. The formation of massive stars and associated 
low-mass population of Trumpler 37 continuing for a few Myr, and the most massive 
star in Trumpler 37, HD 206267, was likely born < 3 Myr ago. Then, the IC 1396 H II 
region was formed by strong UV radiation. In the last stage, the H II region presumably 
triggered the formation of low-mass stars in the dense inhomogeneity. Our results may 
suggest that massive stars were born after the continuous formation of low-mass stars for 
10 Myr. The birth of the exciting star as the most massive one could be the late stage of 
slow but contiguous star formation in the natal molecular cloud. It presumably triggered 
the formation of many low-mass stars at the dense inhomogeneity in and around the H II 
region by a radiation-driven implosion. Spectroscopic followup of our PMS candidates to 
establish their membership correctly and identify their nature is needed to clarify the star 
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formation history of this region. 
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Fig. 1. — Boundary of the observed fields shown on th e Digitized Sky Surve y R image of IC 
1396 . Two O stars, bright- rims of IC 1396 fRi m A-H |Pottaschlll956l : SFO: Isugitani et al. 



1991 



and dark clouds (B 



Barnard 



19271 : Kh 



Khavtassi 



1960) are labeled by their name. 
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Fig. 2. — (a) Spatial distribution of the detected emission-line stars in IC 1396. (b) Contour 
map of the surface density of emission-line stars. The isolines are drawn at intervals of 150 
stars per deg 2 , starting from 100 stars per deg 2 . The position of HD 206267 is shown by a 
star symbol near the center, and the boundary of the observed fields is shown by a solid line. 



-33 - 



S 5 5 

cd 



CO 4 
-Is 3 

2 
1 




CO 

^3 
CO 

-D 2|- 



1 - 



1 


1 

1 
1 


T 1 1 1 

fa) 


575 
..309 


1 


436CHH) 

• 


.0Q7 / 

/ 


• 




( 


• 

• • 

• 

• 
• 

• 


• 


1 


2 

S9W 


3 4 5 6 

- 1 1 RW rmanl 
l i ovv [ii idy j 


r 


r 


1 1 

lb) 




575 

• 

• 


309 436(HH) " 


r 

/ 

/ 287 

1 m 

1 


7 • 

/ 

; 
i 

- : 


i 

1 

1 

1 

1 




i 

* / . 


1 

1 

1 

1 

1 

1 

1 

1 

• 1 

1 


• • 
• 

• 


/ 

• / 

•/ 

; 
/ 
/ 


■ i i i 



Ks - S9W [mag] 



8 10 




2 4 6 8 

Ks - L1 8W [mag] 



10 



Fig. 3.— (a) (J-Ks) vs. (S9W-L18W), (b) (J-Ks) vs. (Ks-S9W), and (c) (J-Ks) vs. 
(Ks-L18W) diagrams of our Ha stars identified in the AKARI/IRC mid-infrared point 
source catalog. Three carbon stars and the candidate for the exciting source of the HH 
588 object are labelled by their corresponding source numbers. The dashed line in (a) is 
the cr iterion separating asymptotic giant branch stars from TTS adopted by Takita et al. 



(2010). The dashed areas in (b) and (c) indicate the areas occupied by TTS (ITakita et al. 
2010). 
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Right ascension 



Fig. 4. — Contour map of the number of emission-line stars overlaid on a false-color map 



using the infrared cam era 



Contreras et al. 



fl2002h and 



RC) on AKARI . Blue is 9 /im and red is 18 /j,m. O-F stars from 



deZeeuwetal 



(119991 ) are also shown. The star symbol near the 
center shows the position of the O star HD 206267. The boundary of the observed fields is 
shown by the white solid line. 
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Fig. 5. — NIR color-color diagram. The thick and thin solid curves are the loci of dwarfs 
and giants, respective ly. The data for the 09-M6 dwarf and G0-M7 gi ants are from 



Bcssc 



1 k Brett 



Meyer et al. 



(119881 ) . The dash-dotted line is the unreddened CTTS locus of 
(119971 ). Emission-line stars are shown by dots, with the large dots representing emission-line 
stars with extinction Av > 1.0. Squares indicate Class-I or Herbig Ae/Be star candidates. 
Three possible carbon stars are shown by open circles. The Av=5 reddening vector is also 



shown. 
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Fig. 6. — Histogram of measured values of source extinction. The vertical line indicates Av 
=1.0. 
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Fig. 7. — Histogram of Ha equivalent width. The shaded inner histogram represents probable 
stellar members of the IC 1396 region. 
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Fig. 8. — Contour map of the number of probable members. The isolines are drawn at 
intervals of 100 stars per deg 2 , starting from 100 stars per deg 2 . The dashed boxes represent 
the subregions referred to in the text. The position of HD 206267 is shown by a star symbol 
near the center, and the boundary of the observed fields is shown by a solid line. 
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Fig. 9. — Extinction-corrected i' vs. {if — J) diagram of probable members. Ha stars 
identified with known T Tauri stars with low extinction are shown by open circles. The 0.1, 



M Q of 



, 10, and 10 


3 Myr 


Siess et al. 


(2000) 



(120001 ) are overlaid. The reddening vector corresponding to j4v=1.0 mag 



is also shown. 
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Fig. 10. — Distribution of ages ( left panels) and ma sses (right panels) of 



from the choice of PMS models (ISiess et al. 
Table 1: Log of Observations. 



(J2000|) and 



i a sta rs predicted 



Palla fc Stahler 



jl99flh ). 



UT Date 


number of fields 






grism 


z'-band 


2004 Nov 12 


1 




2006 Sep 6 


2 




2006 Nov 18, 20 


5 




2007 May 7 


2 




2008 Aug 1, 2, 4, 6, 7 


52 


54 


2008 Aug 8,9,11 


38 




2008 Aug 21, 22, 24, 26 


34 


80 


2008 Nov 2 


2 




2009 Aug 13, 15, 17, 18 


21 


54 


total 


157 


188 a 



a including overlapped areas. 
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Fig. 11. — Distribution of (a) ages and (b) masses of Ha stars in each of our selected areas. 
The ages and masses are derived from a comparison with the Siess model in Figure 9. The 
SF037/38/39/41/Rim B regions are combined with the BRC region. 
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Fig. 12. — Spatial distribution of Ha stars of age (a) < 1 Myr, (b) 1-3 Myr, and (c) > 3 
Myr. Class-I or Herbig Ae/Be stars selected from Figure 5 are also shown by crosses in the 
left panel of (a). Stars with ages 3-5 Myr, 5-10 Myr, and > 10 Myr are shown by open 
circles, open squares, and crosses in the right panel of (c). The star symbol near the center 
shows the position of HD 206267. 



Table 2. List of emission-line stars in IC1396 



^o 




RA 


Dec 




EW[A] 


2MASS b 


quality flag 


J c 


J-H c 


H-Ks c 


Av d 






(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


i 


21: 


:30:53.9 


58:05:20.2 


14.59 


3 


21305396+5805201 


AAA 


13.17 


0.41 


0.33 




2 


21: 


:32:07.1 


57:17:01.7 


15.85 


8 


21320707+5717018 


AAA 


14.13 


0.49 


0.12 




3 


21: 


:32:08.9 


57:44:22.3 


- 


5 


21320901+5744224 


AAA 


14.42 


0.72 


0.34 


0.14 


4 


21: 


:32:12.2 


57:05:40.4 


15.33 


5 


21321218+5705406 


AAA 


14.13 


0.23 


0.19 




5 


21: 


:32:18.4 


56:54:51.8 


16.49 


9 


21321834+5654517 


AAA 


14.65 


0.57 


0.10 




6 


21: 


:32:25.6 


57:49:40.3 


17.02 


46 


21322559+5749403 


AAA 


14.38 


1.08 


0.60 


3.02 


7 


21: 


:32:36.9 


57:57:27.8 


15.30 


3 


21323682+5757276 


AAA 


13.84 


0.54 


0.16 




8 


21: 


:32:49.1 


57:13:48.3 


15.06 


3 


21324900+5713478 


AAA 


13.18 


0.66 


0.15 




9 


21: 


:32:52.4 


57:19:08.5 


18.54 


2 














10 


21: 


:33:01.5 


57:48:34.4 


15.05 


12 


21330146+5748338 


AAA 


12.69 


0.91 


0.29 





(12) 



GO 



Note. - - Table 2 is published in its entirety in the electronic edition of Astronomical Journal. A portion is shown here for 



guidance regarding its form and content. 

a "s" and "-" denote saturated and low photometric quality, respectively. 
b "6x" represents 2MASS 6x Point Star Catalog 
C 2MASS system. 

d calculated from the NIR color-color diagram 

e The "Class" of the object, Class-I (I), Class-II (II), Class-Ill (III) , transition object (TO), Herbig Ae/Be star (HAeBe), 
mainly from Sicilia-Aguilar et al. (2006b) are shown in the parenthesis. VES Coyne et al. (1983), [OH83] Ogura & Hasegawa 
(1983), MVA Marschall & van Altena (1987), [SBZ] Schulz et al. (1997), CGCS Alksnis et al. (2001), [OSP] Ogura et al. 
(2002), [RRY] Reach et al. (2004), [SHB] Sicilia-Aguilar et al. (2004), [SHH] Sicilia-Aguilar et al. (2005), [SHF] Sicilia-Aguilar 
et al. (2006b), [SHW] Sicilia-Aguilar et al. (2007), [GFG] Getman et al. (2007), [ISW] Ikeda et al. (2008), [MMC] Mercer et 
al. (2009), [BVD] the serial number of T Tauri candidates in Table Bl of Barentsen et al. (2011). 

f carbon star candidate 



s single 2MASS source 
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Table 3. AKARI and 2MASS photometric data of Ha stars identified by AKAR1 PSC. 



No 




[L18J¥] a 


J-Ks 


Ks 


remarks' 3 


13 


4.44 


2.88 


2.42 


9.83 


LT Cep, IRAS 21317+5734 


21 


5.72 


3.91 


1.20 


10.64 




29 


4.62 


2.69 


2.21 


9.75 




59 


6.29 


4.24 


1.44 


10.22 


V488 Cep 


60 


5.81 




1.94 


12.49 




81 


6.99 




1.66 


10.54 




115 


5.88 


3.85 


2.53 


9.39 


[RRY] tet 


121 


6.58 


4.50 


0.46 


14.42 




127 


6.77 


4.89 


1.58 


10.03 




203 


5.87 


3.81 


0.73 


9.84 




236 


5.93 


3.79 


1.74 


8.77 


MVA 426 


252 


4.99 


3.04 


1.69 


8.59 


GM Cep 


287 


4.64 


4.06 


2.79 


6.96 


Kun 80, CGCS 5401 c 


309 


2.94 


2.20 


3.56 


6.94 


c 


oou 


u.oy 


A 

^±.00 


n qq 


iu.uo 




377 


6.47 


4.73 




10.76 


[SHF] 93-361 (II) 


436 


4.29 


0.11 


3.50 


12.80 


IRAS 21388+5622, HH 588 


456 


3.80 




2.29 


10.67 


[OSP] BRC38 9 


487 


7.43 




0.72 


8.07 


Kun 332 


575 


1.17 


0.47 


3.68 


4.83 


CGCS 5434 c 


596 


6.80 


5.11 


2.07 


11.24 
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Table 3 — Continued 



No 


[S9W] a 


[L18W] a 


J-Ks 


Ks 


remarks' 3 


606 


7.07 




2.72 


11.67 




611 


6.61 


3.81 


1.64 


9.66 


IRAS 21440+5648 


617 


7.07 




1.49 


10.18 


Kun 215 


618 




5.03 


1.76 


11.32 




621 


6.26 


3.62 


2.02 


10.63 




622 


6.63 


3.68 


1.59 


10.87 





a The zero-magnitude flux densities for the AKARI bands are 
56.262 Jy and 12.001 Jy for the S9W and L18W bands, respec- 
tively. 

b MVA Marschall & van Altena (1987), CGCS Alksnis et al. 
(2001), [OSP] Ogura et al. (2002), [RRY] Reach et al. (2004), 
[SHF] Sicilia-Aguilar et al. (2006b), [GFG] Getman et al. (2007) 

c carbon star candidate 
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Table 4: Summary of the numbers of Ha stars in the surveyed area. 

sky area number of Ha stars a number density 





deg 




per c 


leg 2 


East 


0.11 


89 (34) 


8.1 


X 


10 2 


West (Rim A) 


0.11 


116(56) 


1.1 


X 


10 3 


SF038 


0.13 


51 (25) 


3.9 


X 


10 2 


SF037 


0.10 


15 (7) 


1.6 


X 


10 2 


SF039/41 


0.34 


43 (16) 


1.3 


X 


10 2 


Rim B 


0.09 


23 (12) 


2.8 


X 


10 2 


Halo 


0.99 


167 (38) 


1.7 


X 


10 2 


Field 


2.33 


135 (17) 


5.8 


X 


10 1 


total 


4.2 


639 (205) 









a The number in parentheses represents the probable members. 12 Ha stars with low extinction, HAeBe/Class- 
I candidates (see §4.1), and No. 436 (HH 588 center) are included as probable members. 
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Table 5: Ages an d masses of Ra stars associated with IC 1396. 



Age/Mass 


East 


West 


BRC 


Halo 


Field 


total 


< 1 Myr 


2 


15 


17 


7 


2 


43 


1-3 Myr 


18 


26 


23 


16 


4 


87 


3-5 Myr 


3 


4 


5 


6 





18 


5-10 Myr 


5 


4 


5 


2 


1 


17 


10 Myr < 


2 


2 


2 


5 


5 


16 


0.1-0.2 M 


1 


5 


6 


2 





14 


0.2-0.4 M 


1 








2 


3 


6 


0.4-0.6 M 


5 


9 


9 


9 


1 


33 


0.6-0.8 M 


7 


7 


15 


11 


3 


43 


0.8-1.0 M 


15 


29 


22 


12 


4 


82 


> 1.0 M 


1 


1 








1 


3 


total 


30 


51 


52 


36 


12 


181 



